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We extend our previously developed coarse-graining strategy for linear polymers with a tunable 
number n of effective atoms (blobs) per chain [D'Adamo et al., J. Chem. Phys. 137, 4901 (2012)] to 
polymer systems in thermal crossover between the good-solvent and the 9 regimes. We consider the 
thermal crossover in the region in which tricritical effects can be neglected, i.e., not too close to the 
6 point, for a wide range of chain volume fractions $ = c/c* (c* is the overlap concentration), up to 
$ « 30. Scaling crossover functions for global properties of the solution are obtained by Monte Carlo 
simulations of the Domb- Joyce model with suitably rescaled on-site repulsion. They provide the 
input data to develop a minimal coarse-grained model with four blobs per chain (tetramer model). 
As in the good-solvent case, the coarse-grained model potentials are derived at zero density, thus 
avoiding the inconsistencies related to the use of state-dependent potentials. We find that the 
coarse-grained model reproduces the properties of the underlying, full-monomer system up to some 
reduced density $ which increases when lowering the temperature towards the 6 state. Close to the 
lower-temperature crossover boundary, the tetramer model is accurate at least up to $ ~ 10, while 
near the good-solvent regime reasonably accurate results are obtained up to "1? ~ 2. The density 
region in which the coarse-grained model is predictive can be enlarged by developing coarse-grained 
models with more blobs per chain. We extend the strategy used in the good-solvent case to the 
crossover regime. This requires a proper treatment of the length rescalings as before, but also a 
proper temperature redefinition as the number of blobs is increased. The case n = 10 is investigated 
in detail. We obtain the potentials for such finer-grained model starting from the tetramer ones. 
Comparison with full-monomer results shows that the density region in which accurate predictions 
can be obtained is significantly wider than that corresponding to the tetramer case. 
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I. INTRODUCTION 

The phase diagram of solutions of linear polymers and 
their large-scale structure (i.e., on the scale of the coil 
size) are well understood and can be explained by in- 
voking the physical picture of polymers in implicit sol- 
vent with effective (solvent-mediated) interactions among 
the monomers^"— When the number L of monomers per 
chain (degree of polymerization) is large, such chains ex- 
hibit scale invariance, meaning that, in suitably rescaled 
units of density, temperature, and chain length, solutions 
of chemically distinct polymers and solvents obey uni- 
versal scaling relations, i.e., global properties of the so- 
lution are characterized by universal exponents and scal- 
ing functions. Central properties are the single-chain ra- 
dius of gyration Rg and the equation of state for the 
osmotic pressure 11. For small concentrations, 11 admits 
the expansion^ 



^n 



where /3 = l/{kBT), Bj, are the usual virial coefficients, 
and c ~ N/V is the polymer density {N chains in a 
volume V). The coefficients Bk depend on all chemi- 
cal details. On the other hand, the dimensionless ratios 



A, = BkR 



-3(A;-1) 



{Rg is the zero-density radius of gyra- 



l + B2C + B^c^ + o{e'), 



(1) 



tion) have a universal limit, i.e., they are independent of 
chemical details, for large values of L. 

Scaling functions and exponents are influenced by 
the nature of the interactions among the monomers. 
Three regimes can be distinguishedi^ i) The good-solvent 
regime, usually observed at high enough temperature, 
in which interactions are dominated by the pairwise re- 
pulsion and the single-chain size grows^^ as L'', where^ 
v = 0.587597(7) ~ 3/5 is the Flory exponent. The 
second virial coefficient is positive andiS A2 = 5.500(3) 
for L — ^ 00. ii) The -regime at lower temperature, in 
which the repulsive and attractive contributions in the 
two-body and three-body effective interactions are van- 
ishingly small and of the same order of magnitude. In 
the scaling limit L — > 00 , this regime collapses to a single 
temperature Tg, called 9 temperature. Kt T = Tg chains 
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FIG. 1: Temperature(r)-concentration(c) phase diagram for 
solutions of polymers of finite length L. 



are Gaussian (ideal) with scaling exponent v = 1/2 and 
the second virial combination A2 vanishes. However, 
for linear chains of finite length, as is the case in ex- 
periments and computer simulations, this regime is ob- 
served in a temperature interval around Tg of the order 
of 1/VX. In analogy with the dilute-gas case, one often 
defines the "Boyle" temperature Tb, as the value of T 
at which the second virial coefficient B2 vanishes. Tb 
is generally larger then Tg, converges io Tg as L ^ 00 
[Tb =Tg + 0(i-i/2(inL)-Vii)]^ and is sometimes used 
to estimate Tg numerically^^ iii) At even lower tempera- 
tures the two-body effective interactions become predom- 
inantly attractive (negative second virial coefficient), the 
chains are segregated by the solvent and below a crit- 
ical temperature Tc the solution separates in a chain- 
poor solution where chains are in the collapsed state, 
and a solvent-poor solution where the chain statistics ap- 
proaches the one in the melt state (ideal at large scales). 
Again, in the scaling limit L — >■ c» the critical tempera- 
ture Tc tends to Tg^-± 

For finite values of i, a gradual transition between the 
good- and ^-solvent regimes is observed, i.e., the observ- 
ables and their associated scaling functions interpolate 
smoothly between the two universal behaviors. A sketch 
of the various regimes in the phase diagram of polymers 
of finite length is shown in Figure [T] At low polymer 
concentration, i.e., for <J> ^ 1, when increasing the tem- 
perature above the 9 point, the coil size is observed to 
behave as L'^ with an apparent exponent that increases 
from 1/2 to 3/5. Analogously, the second virial combi- 



nation A2 varies smoothly in the range < A2 < 5.5, 
with the value 5.5 (plus finite L corrections) obtained in 
the fully developed good-solvent regimeJ^ At larger poly- 
mer concentrations and for any T > Tq, a region where 
2-body and 3-body terms are of the same order of mag- 
nitude is present and is generically referred as tricritical 
region. For small concentrations this region extends in a 
small temperature interval^ii whose width is of the order 
of 1/vL (more precisely, it scales as L~^/^(lnL)~^/^^). 
As the concentration c increases, this region widens as 
c/lnL. In the tricritical region scaling corrections de- 
cay very slowly, as inverse powers of logLA^ ^^^^^^ Out- 
side the tricritical region the physics is again dominated 
by two-body effects and by the unbalance between the 
short-range repulsion and the long-range attraction, as it 
was shown by computer simulations of physical polymer 
models.— Here we refer to the thermal crossover region, 
as the region in the phase diagram between the tricrit- 
ical and good-solvent regimes. At finite concentration, 
scaling argumentsi^ indicate that the size of the thermal 



crossover region grows with c as c 



l/(6i/-2) 



^0.40 



At 



sufficiently large temperatures, the good-solvent regime 
is reached at any concentration. 

In the thermal crossover region there are experimental 
and computer-simulation evidences that, for sufficiently 
large values of L, global properties of the solution follow, 
up to small corrections in T, $, and L, general relations 
of the form 



0(T,c,L) = «iOg(^,c)/o(z,$), 



(2) 



where Og{L,c) is the expression of O for the Gaus- 
sian chain, the function fo is called crossover function, 
z = a2(T — Tg)L^'^, and ai and a2 arc nonuniversal 
constants that embody all chemical details. Theory^i^ii^ 
supports the scaling behavior ([2]), albeit with a slightly 
different scaling variable. Indeed, the crossover limit 
should be taken by keeping a2{T - Tg)L^/^{\nL)-^^^^ 
fixed, which differs by a power of In L from the scaling 
variable z = {T — Tg)L-^/^ which appears in Eq. ^. Such 
a logarithmic dependence is irrelevant in all practical ap- 
plications, since the observation of this slowly varying 
term would require data in a very large interval of poly- 
mer lengths/molecular weights. By varying z one obtains 
the full crossover behavior, from the tricritical region, 
corresponding to small values of z, to the good-solvent 
regime, which is obtained for z — > 00. 

It is important to note that the nontrivial universal be- 
havior is obtained when taking simultaneously the limits 
L — ^ 00 and T — >■ Tg in such a way that the arguments 
of the crossover function fa, z = a2{T — Tg)L^''^ and 
$, remain constant. If the limits are taken differently, 
one would obtain a different result. For instance, if one 
takes the limit L — > 00 at fixed T > Tg, one would obtain 
good-solvent behavior in all cases, while, if one decreases 
T towards Tg at fixed large L, only the tricritical be- 
havior would be observed. A similar variety of universal 
scaling behaviors, which depend on how the limit L — )■ 00 
is taken, is found in the stretched-chain problcmJ^ 



The variable z that parametrizes the temperature 
crossover depends on the nonuniversal parameter a^^ 
hence a proper definition requires specifying a physically 
meaningful normalization condition. To avoid this prob- 
lem, one can proceed as suggested in Refs. ll4lll7l . i.e., one 
can parametrize the crossover in terms of a physical vari- 
able rather than in terms of z. For instance, one can use 
the dimensionless sccond-virial combination Ai- Then, 
Eq. ([2]) can be written as 



0(T,L,c) = aiOG(i,c)5o(^2,<i>) 



(3) 



where 50(^21 ^) is universal. The quality of the solution 
is now characterized by A2 that varies between zero at 
the Q point and ^2,05 — 5.5, the good-solvent valueJ^ In 
the thermal crossover region the relation between 2 and 
Ai has been obtained combining analytical results and 
Monte Carlo simulationsr^ obtaining 



Ai{z) ^ 4tt^/'^z{1 + 19.1187Z 
+126.783z2 + 331.99z3 



(4) 



268.96Z'' 



-1/4 



where z has been normalized so that ^2(2;) w 47r^'^z for 
small z, i.e., *(z) « z, where * = 2{4:Tr)~^/'^A2 is the 
so-called interpenetration ratio, which is commonly used 
in much of the experimental literature. 

As discussed by Sokal^ the universal scaling func- 
tions that parametrize the thermal crossover coincide 
with the crossover functions that are defined in the two- 
parameters model (TPM) (ii^ii^ which considers only re- 
pulsive two-body monomer-monomer interactions. In 
this framework the scaling variable z is identified with the 
usual Zimm-Stockmayer-Fixmann variablci^ Therefore, 
to determine the leading crossover behavior, it is con- 
venient to work in the TPM framework, which provides 
directly the crossover functions. This is the approach we 
have taken in this paper, using the lattice Domb- Joyce 
modelfSi which is the lattice version of the continuum 
TPM, to obtain the crossover functions, both in the di- 
lute and in the semidilute regime. 

Although the physics of homogenous polymer solutions 
is well understood, in more complex situations — for in- 
stance, inhomogeneous cases or when the polymers are 
only a single component of a more complex system — it 
might be important to reduce the number of degrees of 
freedom representing the polymer subsystem or, cquiv- 
alcntly, to limit our interest to the physics occurring at 
a length scale comparable to the polymer size. In these 
cases adopting a coarse-grained (CG) representation of 
the solution might be essential for the feasibility of the 
numerical investigation. Moreover, since the target prop- 
erties on which the coarse-grained model (CGM) is built 
can be determined in the scaling limit i — > 00, CGMs 
allow one to study thermodynamics and structural prop- 
erties in the universal, scaling limit without requiring ad- 
ditional extrapolations. 

Several routes to CGMs can be adopted. In the 
structure-based route22ri2^ the CGM is set up in such 



a way to reproduce the marginal probability distribu- 
tion of a set of chosen structural collective variables. 
This procedure, however, does not predict the cor- 
rect thermodynamics of the underlying solution. Con- 
versely, a CGM can be defined to match the ther- 
modynamic behavior, but then structural properties 
are not reproduced correctly»2^ A third method is the 
force-matching approach (often called multiscale coarse- 
graining method^l), in which the state-dependent pair- 
wise potential is determined by requiring the CG sys- 
tem to match the atomistic force on the CG atoms as 
accurately as possible. Also this method has a struc- 
tural interpretation)^ the matching condition is equiv- 
alent to require the CG force to satisfy the appropriate 
Yvon-Born-Green equation that relates the pair and the 
three-body correlation function. However, no guarantee 
of reproducing the pair structure and the thermodynam- 
ics is provided. The inconsistency between structure and 
thermodynamics stems from the fact that CGMs neglect 
multi-body interactions, which would be required to ob- 
tain an exact mapping of the microscopic model onto the 
CGM. 

Recently, we have introduced a procedure to set up 
a hierarchy of CGMs for linear chains in good-solvent 
conditions which simultaneously reproduce quite accu- 
rately structure and thermodynamics of polymer solu- 
tions deep into the semidilute regime i^^'^° The minimal 
model consists in representing each linear chain by a short 
polyatomic molecule with four CG sites (tetramer). The 
tetramcr potentials are set up at zero density — this allows 
us to avoid the inconsistencies^"— that occur when us- 
ing state-dependent potentials — by matching the single- 
chain intramolecular structure and the center-of-mass 
pair correlation function between two identical chains. 
The minimal representation (tetramcr) has been shown 
to provide accurate results for the underlying solutions 
up to $ ~ 2i^ Larger values of $ can be accessed by 
increasing the resolution of the CGM, i.e. by increasing 
the number of effective monomers (blobs) per chain. The 
effective potentials for these higher-resolution models are 
obtained by using a simple transferability approach. We 
assume that the potentials are independent of the number 
of blobs, as long as all lengths are expressed in terms of 
the blob radius of gyration. This transferability approx- 
imation was shown to be quite accurate^ and allowed 
us to obtain accurate thermodynamic and (large-scale) 
structural results for $ ::^ 1. 

In this paper we present the extension of this CG strat- 
egy to solutions in the thermal crossover region, defin- 
ing CGMs with z-depcndcnt potentials that reproduce 
the polymer crossover functions. We have considered the 
TPM as the underlying microscopic model, using the lat- 
tice Domb- Joyce (DJ) model^l as reference system. As 
in the good-solvent case, we first determine the effective 
interactions for the tetramer case for several values of 
z in the crossover region, by matching structural prop- 
erties at zero density. These potentials are then used 
to define CGMs of higher resolutions, i.e., higher num- 



ber of CG sites. The transferability approach is more 
complex than in the good-solvent case, since one must 
change at the same time the reference length scale and 
the scaling variable z. We perform extensive simulations 
of the CGMs at finite density, both in the dilute and 
in the semidilute region. The results are then compared 
with those obtained in full-monomer (FM) simulations 
and with the field-theory predictions of Rcf. y. We find 
that the tetramer CGM reproduces very accurately the 
crossover functions up to $ « 30 (and quite reasonably 
at all densities) for small values of z. Close to the good- 
solvent behavior, instead accurate results are obtained 
up to <& ~ 2. Use of the decamer model with n = 10 
blobs allows one to widen significantly the density region 
in which the CGM is predictive. 

The present work represents an important extension 
of our previous results, allowing us to consider polymer 
systems in the thermal crossover regime. In particular, 
it opens the way to the study of more complex systems 
like polymer-colloid solutions away from the good-solvent 
regime and diblock copolymer solutions,—!^ where each 
block has different affinity with the solvent. 

The paper is organized as follows. In section |ll] we 
define the DJ model we have adopted and explain how 
the crossover functions are computed. In section IIIII we 
present our CG strategy, first illustrating the CG repre- 
sentation (CGR) of the FM chain fsection FlII A[) and later 
presenting our CGM (section IIIIB[) . Section IIVI is de- 
voted to the comparison between the results of the CGR 
of the FM system and of the CGM, first at zero density 
(section IIV A[) and then at finite density in the semidi- 
lute regime (section IIVC[) . Simulation results are also 
compared with field-theoretical expressions^ and large- 
density predictions obtained by using the random-phase 
approximation (RPA), which becomes exact for $ — )■ oo. 
In section |Vl the transferability of the tetramer poten- 
tials to finer resolution models is illustrated and validated 
against FM predictions. Finally, section IVTl reports our 
conclusions. Two appendices are also present. In Ap- 
pendix [X] we report the calculation of the blob radius of 
gyration in the crossover regime, while in Appendix [Bj we 
summarize the field-theory predictions of Ref. |5|- Some 
numerical details and explicit expressions for the CGM 
potentials can be found in the supplementary materiali^ 



II. CROSSOVER FUNCTIONS FROM MONTE 

CARLO SIMULATIONS OF THE DOMB-JOYCE 

MODEL 

In order to compute the TPM crossover functions, 
we consider the three-dimensional lattice Domb- Joyce 
(DJ) modeli^ In this model the polymer solution is 
mapped onto N chains of L monomers each on a cubic 
lattice of linear size M with periodic boundary condi- 
tions. Each polymer chain is modelled by a random walk 



nian is given by 



JS) 



JS) 



{ri^...,r^'^}with|r^'^-r 



(») JS) 



1 (we take the lattice 



N 



i=l l<a</3<L 

+ E EE^(ri^^4^')' 



(5) 



where (5(r, s) is the Kronecker delta. Each configuration 
is weighted by e^™^, where w > is a free parameter 
that plays the role of inverse temperature. This model is 
similar to the standard lattice self-avoiding walk (SAW) 
model)^'^ which is obtained in the limit w — > +oo. For 
finite positive w intersections arc possible although ener- 
getically penalized. For any positive w, this model has 
the same scaling limit as the SAW model^i and thus al- 
lows us to compute the universal scaling functions that 
are relevant for polymer solutions under good-solvent 
conditions. 

The TPM results in the thermal crossover region can 
also be derived from simulations of the DJ model. They 
are obtained^ by taking the limit lu — > 0, L — > oo at 
fixed X = wL^I^. The variable x interpolates between 
the ideal-chain limit {x = 0) and the good-solvent limit 
[x = oo). Indeed, for w = the DJ model is simply 
the random- walk model, while for any w 7^ and L — > 
00 one always obtains the good-solvent scaling behavior. 
The variable x is proportional to the variable z that is 
used in the TPM context. If we normalize z so that 
^2(2) ~ A-TT'l'^z for small z as in Eq. (H]), we have^2ii£ 



3/2 



w 



L^'\ 



(6) 



As discussed in Ref. Il5|, the TPM results can be obtained 
from Monte Carlo simulations of the DJ model by prop- 
erly extrapolating the numerical results to L — >■ 00. For 
each z we consider several chain lengths Li. For each 
of them we determine the interaction parameter Wi by 
using Eq. ©, that is we set Wi = (27r/3)^/^zL~"^'^. Sim- 
ulations of chains of Li monomers are then performed 
setting w = Wi. Simulation results are then extrapolated 
to L — > C!0, taking into account that corrections arc of 
order 1/vL.— 1^ If R{L, z) is a dimensionless ratio of 
two global quantities at zero density (for instance, the 
second- virial combination A2), the TPM result R*{z) is 
obtained by performing an extrapolation of the form 



R{L, z) = R*{z) + ^+ 0{L-^ InL 

vL 



(7) 



spacing as unit of length) and 1 < i < N. The Hamilto- 



Of course, for z — >■ 00, R*{z) converges to its universal 
good-solvent value, which can be obtained by taking the 
limit L — >■ 00 at fixed (arbitrary) w. 

In this work we often consider adimensional distribu- 
tion functions g(p] L, z) that also depend on the adimen- 
sional ratio p = r/Rg. For these properties, the TPM 



TABLE I: Values of z considered in this paper. Tlie notation 
is such that ^2(2'^') « p. The estimates of yl2(z) for z^^^ 2*^^ 
z^ \ z'*', and z'^', are the direct MC estimates of Ref . [iSl: in 
the other cases we use the interpolation formula Q. In the 
last column we report r = Aiiz^jA^^GS ~ ^(•2)/^GS {'^ is 
the interpenetration ratio often used in experimental work), 
where A2,gs and '^gs are the good-solvent values. 



Mz) 



.(1) 



j(1.5) 
.(2) 



0.056215 0.9926(10) 0.18 
1.5 0.27 

Z-' = 0.148726 1.9782(18) 0.36 

^(2.5) ^ 0.225292 2.5 0.45 

z^^^ = 0.321650 2.9621(27) 0.54 

^(3.5) ^ 0.493088 3.5 0.64 

z'"*) = 0.728877 3.9433(34) 0.72 

^(4.5) ^ 1,32527 4.5 0.82 

z'-^'> = 2.50828 4.9147(36) 0.89 



result g*{p;z) is obtained by performing an extrapola- 
tion of the form 



g{p;L,z) ^ g*{p;z) + 



Vl 



0{L-^\nL). (8) 



At finite density we should take the limit i — > 00, keeping 
the polymer volume fraetion $ fixed. In practice we keep 
the dimensionless combination $ = ATTc[Rg{Li,Wi)]'^ /3 
fixed, where Rg{L,w) is the zero-density radius of gy- 
ration (here and in the following we indicate any zero- 
density quantity with a hat). 

In order to determine the crossover behavior, we have 
performed simulations at the five values z^*\ i = 1, . . . 5, 
considered in Ref. [l^, see Table HI They belong to the 
crossover region between ideal and good-solvent behavior 
and are such that ^2(2^")) « n (remember that ^2(2) 
varies between and 5.50). Some additional simulations 
have also been performed at four values of z such that 
A2{z) = 1.5, ..., 4.5, where ^2(2) is given in Eq. (g]). 
The explicit values are reported in Table HI 

In principle, the crossover functions can also be com- 
puted in other models that interpolate between good- 
solvent and 6 behavior. A typical example, which 
has been widely discussed in the literature, see, e.g., 
Refs. [IIII41II42I is the SAW model with nearest-neighbor 
interactions (often called interacting SAW model). An- 
other model, quite interesting from a computational 
point of view, is the extension of the D J model discussed 
in Ref. Il3l in which an additional energy term associ- 
ated with triple intersections is considered. Unlike the 
DJ model, in these models there is a nonvanishing three- 
body effective coupling, which implies, for instance, that 
the third virial combination A3 is positive at the Boyle 
temperature where A2 = 0. If we take these marginally 
irrelevant three-body terms into account, the scaling be- 
havior ^ becomes^i^ii^ 



where U3, which parametrizes the effective three-body 
interaction, vanishes as 1/lnL for L — ^ 00. As long as 
the thermal crossover region is considered, i.e., we are 
in the temperature/degree-of-polymerization region such 
that z 3> U3 and L is large, we can neglect uy, in Eq. (|9]) 
and rcobtain Eq. ^ with fo{z,<^) = /ic)(z,0,$). How- 
ever, for polymers of finite length, tricritical corrections 
proportional to W3 give rise to slowly varying scaling cor- 
rections, which make a precise determination of the scal- 
ing crossover functions quite difficult. For this reason it 
is convenient to consider models, like the DJ one, with 
only two-body repulsion, thereby avoiding unwanted tri- 
critical corrections. Of course, if one wishes to discuss 
also tricritical effects, the model of Ref. [ij or interact- 
ing SAWs should be considered. In this respect, we note 
that the extended DJ modelH is computationally much 
more convenient than the more common interacting SAW 
model. Indeed, since interactions are soft, one expects 
the Monte Carlo dynamics to be significantly faster than 
for interacting SAWs. 



III. THE BLOB MODEL 

A. The coarse-grained representation of the 
polymer model 

In the multiblob approach one starts from a coarse- 
grained representation (CGR) of the underlying full- 
monomer (FM) model, which is obtained by mapping 
a chain of L monomers onto a chain of n blobs, each 
of them located at the center of mass of a subchain 
of m = L/n monomers. If the monomer positions are 
given by {ri, . . . , tl}, one first defines the blob positions 
Si, . . . , s„ as the centers of mass of the subchains of length 
TO, i.e. 



S,; = 



^ mi 

- E 



(10) 



Q— 7n(z— 1)-|-1 



For the new CG chain {si,...,s„} one defines several 
standard quantities. First, one defines its square radius 
of gyration 



R 



g,b 



n) = — TT 



1 
2~2 



jj=i 



(S,; Sj f 



(11) 



Such a quantity is always smaller than R^ , because of the 
exact identity 



Rl^Rlb{n)+rl{n), 



(12) 



0{T,L,c) = aiOG(L,c)/io(z,U3,$), 



(9) 



where rg (n) is the average radius of gyration of the blobs. 
The ratios R^ ,{n)/R^ and r'j^{n)/R^ of their averages 
over the polymer configurations (we use the same sym- 
bol both for the radius of a single chain and for its aver- 
age; the correct interpretation should be clear from the 



context) show a universal erossover behavior, i.e. inde- 
pendent of the nature of the underlying polymer model 
as long as L and m are large enough. This crossover 
can be equivalently parametrized in terms of z or of A2. 
Explicit zero-density results are reported in App. [X] 

To define the CGM we proceed as in our previous 
work )^^i'^° defining at first a tetramer model with n = 4 
blobs. Higher-resolution models with n > A will be dis- 
cussed in Sec. El To determine the four-blob CGM, we 
compute several intramolecular CGR structural distribu- 
tions for an isolated chain (zero-density limit). First, we 
determine the bond-length distributions of the CGR FM 
model: 



TABLE II: Estimates of R^jj/R^ and of A2 obtained by using 
the tetramer (t) CGM and the full-monomer (FM) model. 
The values of 2 are reported in Table |T] 





Rq,blRg 


A2 = B2/RI 


z 


FM t 


FM t 



,(1) 



,(3) 
^(4) 



.(5) 



0.7553(7) 0.7574(3) 0.9926(10) 0.9763(1) 

0.761(1) 0.7587(5) 1.9782(18) 1.9459(3) 

0.7686(15) 0.7719(5) 2.9621(27) 2.9364(5) 

0.777(2) 0.7844(4) 3.9433(34) 3.9999(7) 

0.787(2) 0.7885(4) 4.9147(36) 4.9105(8) 



P,,(r) = (5(|s,-s, I -r)). 



(13) 



In the crossover limit at fixed z, the adimensional combi- 
nation RgPij{r) converges to a universal crossover func- 
tion f*Ap]z), p = r/Rg, which can be computed in the 
DJ model as described in Sec. HT] Second, we will need 
the distributions of the two equivalent bending angles /3i , 
of the torsion angle 6*, and of the angle /3i3, defined in the 
CGR of the polymer model. They are defined as follows: 



cos Pi 



cos Pi^i 



cos 9 — 



h, b 



i+l 



|b,||bi+i| 
bi ba 



Ibillbal' 
(bi X ba) ■ (ba x bg) 
Ibi X ballba x baj 



(14) 
(15) 
(16) 



with bi = Si_|_i — Si. In the crossover limit they converge 
to the TPM distributions f^{cosl3; z), f^ -^^^{cos (3i3;z), 
and/*(0;z). 

Finally, to determine the intermolecular CGM interac- 
tions we will make use of the center-of-mass intermolec- 
ular distribution function. It is defined by 



gcAiir) = (e-'^^-)o,r 



(17) 



where (•)o,r indicates the average over two isolated poly- 
mers, the centers of mass of which are in the origin and in 
r, respectively, and U12 is the intermolecular energy. In 
the crossover limit L — >■ cxj at fixed z, gcM{i') converges 
to a universal function g^jj^jip] z). 



B. The coarse-grained blob model 

The CGM consists of polyatomic molecules of n atoms 
located in {ti, . . . ,t„}. All length scales arc expressed 
in terms of the zero-density radius of gyration, hence 
all potentials and distribution functions depend on the 
adimensional combinations p = t/Rg. In order to have 
an exact mapping of the CGR of the polymeric system 
onto the n-blob CGM, one should consider an n-body 
intramolecular potential, which, for n > 2, can be ex- 
pressed in terms of 3{n — 2) scalar combinations of the 



positions of the blobs because of rotational and trans- 
lational invariance. Even for n as small as 4, this re- 
quires considering a function of 6 independent variables. 
Moreover, since we are computing the potentials in the 
crossover region, we should additionally consider them 
as a function of the crossover variable z. Of course, this 
is far too complex in practice. Hence, we have used two 
different simplifications i^i^ First, we use a limited set of 
interactions. The intramolecular interactions have been 
modeled by introducing seven different potentials, each 
of them depending on a single scalar variable. This choice 
is arbitrary, but, as we have already verified in the good- 
solvent case^Si^ it is particularly convenient and works 
quite well. Second, we have computed the potentials only 
for five different values of z, z^^\ . . ., z^^\ reported ex- 
plicitly in Table ID For other values of z we use a simple 
interpolation, which, as we shall discuss, works quite pre- 
cisely. 

As in Refs. [29li30l . we consider a set of bonding pair 
potentials: blobs i and j of the tetramer interact with 
a pair potential Vij{p;z) with p = |tj — tj\/Rg. Be- 
cause of symmetry we have Vi^^p; z) ~ V2a[p',z) and 
Vi2{p;z) — V34(p;z), so that there are only four inde- 
pendent potentials to be determined. Then, we con- 
sider a bending-angle potential Vh(cos/3;z), a potential 
V6,i3(cos/3i3; z) and a torsion-angle potential Vt{9;z), 
where /3, /?i3, and 9 are defined as in Sec. lIII Al 

To determine the CGM intramolecular potentials, we 
require the CGM to reproduce the adimensional TPM 
bond distributions / *,-(r; z ) and the TPM angle distribu- 
tions defined in Sec. IIII Al To obtain these distributions, 
we perform simulations of the DJ model with chains 
of lengths L = 1000, 2500, and 5000 monomers, using 
the pivot algorithmii^Ti^ Then, we extrapolate the rel- 
evant quantities by using Eqs. (O and ([8]). Once the 
TPM distributions are known, the potentials are ob- 
tained by applying the Iterative Boltzmann Inversion 
(IBI) scheme,22i23i9 

The method works quite precisely: The CGM repro- 
duces quite well the target distribution functions. For 
instance, in Fig. [2] we report the intramolecular distribu- 
tion function (/intra (p) computed by using the CGR of the 
polymer model and the tetramer model for the five values 
of z we consider. The agreement is excellent. Excellent 
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FIG. 2: Intramolecular distribution function p^ gintmip', z) . 
We report tetramer (t) and full-monomer (FM) results at 
zero density for the values of z reported in Table [H For sake 
of clarity, results at different values of z are shifted upward 
according to the legend. 



sion potential becomes very small for small values of 
z: the tetramer becom.es more flexible as the excluded- 
volume interactions become less effective. For instance 
^[Vt(7r;z(i)) - 14(0; z(i)] « 0.01. The four-body bend- 
ing potential Vf,_i3(cos/3i3; z) changes only in the re- 
gion in which cos/3i3 > 0. In particular, by increas- 
ing z, Afc^a = /3[Vb^i3{l; z) - Vb.i3{0;z)] decreases, sig- 
nalling that more elongated conformations arc preferred 
as excluded- volume effects increase. 

As for the intermolccular potentials, we have made 
again some drastic simplifications. First, in the spirit of 
the multiblob approach, we neglect interactions among 
three or more tetramers. The potential between two 
tetramers is still a function of the 6(n — 1) relative po- 
sitions of the blobs, which is again far too complex. We 
have thus simplified the model by considering a single 
intermolccular central pair potential W{p]z): all blobs 
interact with the same potential, irrespective of their 
positions along the tetramer. As shown in Refs. [29ll3(J 
this drastic simplification works quite well. Such a po- 
tential has been obtained by requiring the CGM to re- 
produce the TPM center-of-mass intcrmolecular distribu- 
tion function g^j^j(p; z). The potential W{p] z) has been 
parametrized as 



l3W{p-z) = ci(z)cxp(-C2(2;)p^), 



(18) 



agreement is also observed for the angle distributions (see 
supplementary material^) . As a second check of the ac- 
curacy of the inversion procedure, we compute the ratio 
Rgb/Rg both in the polymer model and in the CGM. 
Again, the results reported in Table |ll] show an excellent 
agreement. 

The potentials are shown in Figs. [3] and HI Ex- 
plicit parametrizations are reported in the supplemen- 
tary material.— The bonding pair potentials Vi2{p]z) 
and V2z{p] z) change significantly with z. First, the po- 
sition pmin of the minimum of the potentials decreases 
from the good-solvent valued /Omin — 0.5 to 0.3 for 
z = z^'^' and is approximately zero for z^^'. This is 
consistent with what one finds^ for the random-walk 
case (z = 0), for which the bonding pair potential is 
a parabola, hence pmin = 0. Second, the potential be- 
comes softer at the origin as z decreases. For instance, 
^12(2) = /3[Vl2(0;z) - Vi2{piain; z)] dccrcascs from the 
good-solvent value 0.8 to as z decreases. More pre- 
cisely, Ai2(z) w 0,0.4,0.7 for z = z''^\ z^^\ z^^\ respec- 
tively. As before, this is due to the fact that excluded- 
volume effects become less relevant as z decreases. Po- 
tential 1^13 (p;z) shows qualitatively the same behavior 
as a function of z. Potential Vu{p]z) has an approxi- 
mate Gaussian shape for all values of z. As expected 
Vl4(0;z) decreases significantly with decreasing z. Since 
excluded- volume effects decrease in this limit, overlaps 
are less penalized. 

In Fig. |4]we report the angular potentials. The bend- 
ing potential changes significantly with z. For large 
z folded configurations with cos/3 ~ 1 are penalized, 
while for small z, they are slightly favored. The tor- 



in terms of two unknown z-dependcnt parameters ci(z) 
and C2(z). They have been determined following the ap- 
proach of Ref. [5l|. Parametrization p8p looks adequate. 
The CGM reproduces quite precisely the FM result, see 
Fig. [S] As a check of the accuracy of the procedure, we 
also compare the second- virial combination ^2- Again, 
the results reported in Table HIl show good agreement (dif- 
ferences are less than 2%). The intermolccular potential 
is shown in Fig. U) As expected, in the crossover region 
from good-solvent to 9 behavior one observes a decrease 
of its strength and a slight decrease of its spatial range. 
The results reported above give the tetramer potentials 
for five different values of z, each of which corresponds to 
a different value A2 of the second- virial combination ^2, 
see Table [J To define the model for all values of z we use 
a simple linear interpolation formula. Given z, we first 
determine A2{z) by using Eq. ([4]) and then an integer 

1 < ni < 5 such that ^"'^ < Mz) < 4"'^^^ (" = 
6 corresponds to the good-solvent case, hence A2 = 
5.500). Then, for each potential we set (cr may be p or 
an angular variable) 



Via;z) 



Mz) 



A 



(m+i) 



2 -y(a;z("i+i)) 



A 



(«i) 



A2{z)-A 



(m+i) 



i("i+i) 



-A 



("l) 



V{a;z^''^>). (19) 



As a check of the accuracy of this interpolation wc have 
chosen four values of z such that A2 — 1.5, 2.5, 3.5, 
4.5, see Table HI For these values we have computed the 




r/Rg(z) 



FIG. 3: Intramolecular tetramer bond potentials for five values of z reported in Table H] points correspond to the numerical 
estimates, while the solid lines are the interpolations reported in the supplementary material.— For the sake of clarity, results 
at different values of z are shifted upward according to the legend. 



tetramer potentials using the interpolation formula (|19|) . 
Then, we have computed again A2 for each value of z 
by using the tetramer CGM. Results are consistent — 
differences are at most 1% — confirming the accuracy of 
the interpolation we use (see supplementary material^ 
for more details). 



IV. COMPARING TETRAMER AND 
FULL-MONOMER PREDICTIONS 



in which each polymer is represented by a monoatomic 
molecule located in the polymer center of mass. The 
SB potentials in the 6'-to^ood-solvent crossover regime 
were computed in Refs. [l4ll4lll43 using the self-avoiding 
walk model. Here we perform the same computation 
more carefully, by using the DJ model. An explicit 
parametrization that satisfies all theoretical constraints 
and reproduces the good-solvent results of Ref. [Ij is re- 
ported in the supplementary materiali^ 



In this section, we compare the CGM structural pre- 
dictions with those obtained for the CGR of the polymer 
model. In Sec. lIV Al wc extend the discussion of Sec HIlBI 
at zero density, while in Sec. lIVCl we consider the semidi- 
lute regime. Here we compare the CGM results with 
the results of FM simulations at finite density, discussed 
in Sec. IIVBI Beside considering the tetramer model, we 
shall also discuss the simpler single-blob (SB) model. 



A. Zero-density results 

By construction, the tetramer CGM reproduces all 
CGR bond-length distributions. However, it is far from 
obvious that the distributions of other structural in- 
tramolecular quantities are correctly reproduced. In 
Fig. ini we report the distribution of the ratio R ~ 
Rg^]j/{H^)^/'^, where Rg_b is the radius associated with 




FIG. 4: Tetramer angular and intermolecular potentials for five values of z reported in Table H] points correspond to the 
numerical estimates, while the solid lines are the interpolations reported in the supplementary material.^ For the sake of 
clarity, results at different values of z are shifted upward according to the legend. 



a single chain, while (i?g) is the average squared radius 
of gyration. The agreement is excellent for all values of 
z wc consider. This is a nontrivial check of the quality of 
the multiblob representation since this distribution is not 
directly related to the bond-length distributions nor to 
those of the intramolecular angles we consider. Clearly, 
the tetramer CGM correctly models shape and size of the 
polymer in the whole crossover region. 

It is also of interest to compare the predictions for the 
virial coefficients defined in Eq. dlj. Since wc matched 
the pair distribution function gcM {'>'), the combination 
A2 = B2Rg^ should be the same in the CGM and 
in the polymer model, as we have already verified in 
Sec. IIIIBI Here we compare the third- virial combination 
^3 — B^Rg^. Differences between the CGM and the 
polymer results allow us to quantify the importance of 
the three-polymer interactions that are not fully taken 
into account by the CGM. Beside tetramer and FM re- 
sults, we also report estimates obtained by using the SB 
model. The results are reported in Table IIIIl The SB 



results underestimate quite significantly the FM results: 
the relative deviations increase from 19% for z — z^^^ to 
30% for z = z'-^K The tetramer results arc significantly 
better. Close to the good-solvent regime (z > z'^^) devi- 
ations arc approximately 3%. As z is lowered, differences 
increase: for z = z^^\ we find a 9% difference. It is im- 
portant to note that, although in all cases the observed 
differences for A3 increase as the 9 point is approached, 
this does not imply that the difference between the CGM 
and FM pressure increases when lowering z. Indeed, per- 



turbation theory gives An/A2 



for z — )• 0. There- 



fore, the third and higher-order virial contributions to 
the pressure become increasingly less relevant as z de- 
creases. As a consequence, as we discuss in Sec. IIVCI 
the CGM provides increasingly more accurate estimates 
of the pressure as the 9 point is approached. 
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TABLE III: Estimates of A>, for the full-monomer (FM) 
model, for the tetramer (t), and single-blob (SB) CGMs. We 
consider values z'"' defined in Table [T] The full-monomer 



.(1) 



,(2) ^(3) 



,(4) 



and 



.(5) 



are the di- 



(FM) results for 

rect MC estimates of Ref. Il5l : in the other cases we use the 

interpolation formula A^i^z) reported in Table III of Ref [la. 



z 


FM 


t SB 


z(i) 


0.0849(8) 


0.077(1) 0.059 


^(1-5) 


0.276 


0.261(6) 0.197 


^(2) 


0.6061(30) 


0.553(8) 0.442 


^(2.5) 


1.15 


1.10(2) 0.867 


^(3) 


1.843(8) 


1.74(2) 1.404 


^(3.5, 


2.89 


2.85(3) 2.264 


^(4) 


4.021(13) 


3.95(4) 3.164 


^(4.5) 


5.72 


5.67(5) 4.569 


^(5) 


7.243(22) 


7.02(6) 5.827 



FIG. 5: Radial distribution functions for a pair of isolated 
chains as a function of the center-of-mass separation p — 
r/Rg. We report tetramer (t) and full- monomer (FM) results 
at zero density for five values of z reported in Table [T] 
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simulations for several values of $ and several values 
of L ranging between 100 and 2000. As in our previ- 
ous work)^ we use a combination of pivot)^"— cut-and- 
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FIG. 6: Distribution of the ratio 7? = Rg^t,/ J {R^). We report 
tetramer (t) and full-monomer (FM) results at zero density 
for the values of z reported in Table [H For the sake of clarity, 
results at different values of z are shifted upward according 
to the legend. 



permute 



52-54 



and reptation moves. Since the penalty 



parameter w is quite small, the algorithm is quite effi- 
cient and for 2*^^^ . . ., z*-*^ we are able to simulate poly- 
mer systems for densities up to $ = 30. For z = z^^\ 
i.e. close to the good-solvent regime, the largest density 
considered corresponds to $ = 20. In the finite-density 
simulation we should also fix the volume V = M^ of the 
box. As in Ref. \5^ we fix the box size so that the number 
N of polymers in the box is never smaller than 100 (typ- 
ically 150 < AT < 1000). We measure the relevant CGR 
intramolecular and intermolecular distribution functions, 
needed for a detailed comparison with the CGM results. 
Moreover, to compare the thermodynamics we compute 



K = 



dc 



(20) 



where c = N/V. This quantity is determined as in 
Ref. 152- Using the compressibility rule^ we can relate K 
to the total structure factor that can be easily measured 
in simulations. If we define 



5(k) 



L^N 



N L 
j=l a=l 



(21) 



then we have 



B. Intermezzo: finite-density tiiermodynaniics by 
full-monomer simulations 



In order to discuss the behavior of the CGM in the 
semidilute regime, we need to obtain FM reference data. 
For this reason we have performed an extensive study of 
the DJ model at finite density. For the five values of z, 
z^^' , ..., z^^' , reported in Table |T] we have performed 



- = lim S{k). 



(22) 



Since we work in a finite box of volume M^, we must 
quantify the finite- volume effects. For fluids they have 
been extensively discussed, see, e.g., Refs. [Sa and refer- 
ences therein. In general, at fixed $ finite-volume quan- 
tities converge to their infinite- volume counterpart (if we 
consider distribution functions or the structure factor we 
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should, of course, take the limit at fixed r or k) with 
corrections of order M~^. We find that these corrections 
are smaller than or, at most, of the same order as the 
statistical errors. To give an example, let us discuss in 
more detail the structure factor S'(k). If we assume 



5'(k;M) =S'(k;oo) 



M3 



OiM- 



(23) 



then, we have 



1 



5(k; M)~S{k- oo) = ^ [S(k; M/2) - 5(k; M)]+0{M-^). 

(24) 
Hence, a rough estimate of the size effects can be ob- 
tained by computing the same quantity for two boxes of 
linear size M and M/2, respectively. To give an idea of 
the effect, let us consider two cases. In the good-solvent 
regime, for $ — 10, L = 250, we have performed sim- 
ulations with M = 32 (AT = 100) and M = 64 (TV = 
802). For k — 2tt/32, the smallest momentum which is 
present in both cases, we obtain S'(k; M) = 0.01319(4) 
and 0.01316(2) for M = 32,64, respectively. Clearly, 
the systematic error on the estimate obtained by using 
M = 64 (the only one we use) is negligible. In the op- 
posite limit, consider z — zi, ^ = 15, L = 2000. For 
k = 27r/64, we obtain 5(k;M) == 0.09831(6), 0.09853(8) 
for M = 64, 128, respectively. Again, Eq. ([M| indicates 
that finite- volume corrections on the largest-lattice result 
are significantly smaller than statistical errors. Hence, to 
obtain the TPM results it is enough to extrapolate the 
finite-i data at the same value of <I> by using Eqs. ([7]) or 

®- . . . 

In the simulation we estimate S'(k) for two different 

wave vectors: ki = (27r/Af,0,0) and ka = (47r/M,0,0), 

where M is the linear size of the cubic box. Then, we 

define^ 



K,: 



"^2 



kl 



kjS{ki) - klS{k2) 



(25) 



where k = 2sin(fc/2), fci = 2Tr/M, and fc2 = An/M. As 
discussed in Ref. [52, the estimator Kest of the inverse 
compressibility converges to K with corrections of order 
M-4, i.e. 



K,st{L,M) = K{L) 



M4 '■ 



(26) 



where we have only included the leading correction. It 
turns out that this term is not negligible, hence an addi- 
tional extrapolation is needed to obtain K{L). Of course, 
Kcst also depends on $ and z, but since they do not 
play any role here — we work at fixed density and z — 
they are omitted. To go further, we must specify the L 
dependence of a{L). Since a{L) corresponds dimension- 
ally to a fourth power of a length, we expect on gen- 
eral grounds and verify numerically in a few cases that 
a{L) ^ Rg ^ L^, with corrections that decay as L^^^^: 
a{L) — L^{b + cL^^/^) . Since the infinite-volume K{L) 
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FIG. 7: TPM estimates of K for different values of z, see Ta- 
ble [T] We report full-monomer (points, FM) and field-theory 
(lines, FT) results. 



should behave as Fi{L, z) in Eq. ([7]), we end up with the 
asymptotic expansion 



A'e.t(i,A/) = A-* 



k L2 / ^ 



(27) 



where K* = K*{z,^) is the TPM resuh. For each $ 
and z we have therefore fitted our simulation results to 
Eq. ()27p . keeping K* , k, b, and c as free parameters, in 
order to obtain the TPM prediction K*{z, $). 

It is interesting to compare the numerical results with 
the field-theoretical (FT) predictions of Ref. (the rele- 
vant formulae are summarized in App.|B]). In all cases we 
observe very good agreement, see Fig. [71 For z = z^^\ 
z(^\ and z^'^', in the quite large range of densities we 
consider, $ < 30, the relative differences are always less 
than 1%. For the two largest values of z, differences are 
somewhat larger but remain below 4%. Note that similar 
conclusions were reached in Ref. |52| for the good-solvent 
case. Hence, the parametrization of Ref. |5| of the FT re- 
sults appears to be quite accurate in the whole crossover 
regime. 

For $ 
to K, 

where ctg^z) is the swelling factor of the zero- 
density radius of gyration, computed explicitly in 



oo and any finite z, K{z,^) converge; 

-3 



15 



(z,^) = kpAiiz)^ with kpAiiz) = 6y/Trzag^{z), 



Ref. 



For 



.(1) 



z^^' we obtain kprniz) 



0.543,1.27,2.32,4.07,8.40, respectively. For each value 
of z, we can compare our result for the largest value of 
$ with this asymptotic prediction. For $ = 30, we ob- 
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tain from simulation K{z,<^) = 16.73(6), 37.1(2), and 
64.8(3), to be compared with Kg^si^, $) = 16.3, 38.1, and 
69.6, for z ~ z*^^', z^'^^ and z'^^\ respectively. Clearly, 
for such large value of <I> the asymptotic formula holds 
approximately only for z < z"'. For larger values of z, 
the linear behavior sets in for values of $ that are signif- 
icantly larger than 30. 

Finally, we should mention that Ref. I15| also gave a 
prediction for Z{z, $) [Eq. (4.24) of Ref. Il5| vahd for 



CO and $ < $n 



c(^) 



,2.53 



Requiring an error of 



at most 5%, the asymptotic expansion was found to be 
predictive for 2 > 2, $ < {z/2)^-^^. The first condition 
implies that we can only consider z = z^^\ while the sec- 
ond condition implies that $ should be relatively small, 
<f> < 2. In this density range, by using the equation of 
state reported in the supplementary material,— we find 
that the asymptotic formula reproduces Z{z^^\^) with 
an error of at most 4%, confirming the correctness of the 
numerical estimates of Ref. Iisl 



C. Comparison in the semidilute regime 

Let us now compare the CGM predictions with the 
polymer CGR results at finite density. Let us begin by 
comparing the intramolecular distribution function. Re- 
sults for z = z^^") and z = z^^^ are reported in Fig. [S] 
In all cases the agreement is excellent, even for $ = 4. 
Similar conclusions are reached for the finite-density dis- 
tribution of Rg^b, see Fig. [5] Clearly, the tetramer model 
correctly reproduces the large-scale structure of the poly- 
mer even in the presence of significant polymer-polymer 
overlap. Note that this is not the case in the good-solvent 
regime 1^ In that case, for $ « 4 the tetramer is more 
swollen than the polymer: the probability for two blobs 
to be at a given distance p < 1 is significantly smaller in 
the tetramer than in the polymer. Analogously, Rgj, is 
typically larger for the tetramer than for the polymer. 

Let us now compare the intermolecular structure. In 
Fig. [TO] we report the center-of-mass distribution function 
9c Alii") for the tetramer, the SB model, and from FM 
simulations. For z = z^^^ both the tetramer and the SB 
model correctly reproduce the polymer structure, even 
for $ = 4. Clearly, many-body interactions are weak, 
hence do not influence significantly the polymer behav- 
ior, even in the presence of significant overlap among the 
polymers. For z = z^'^' the SB model reproduces gcAiir) 
for $ = 1, while discrepancies are observed for $ = 4. 
On the other hand, the tetramer model is quite accurate, 
even for $ = 4. Many-body interactions are relevant, 
but are reasonably taken into account by the tetramer 
model. 

Let us finally, discuss the thermodynamics, comparing 
the estimates oi K{z, $) for polymers with those obtained 
by using the SB and the tetramer model, see Fig. [TT] 
(see supplementary material^^ for tables of data). For 
z = z*-^^ tetramer results are on top of the FM results, 
even for $ = 10: at this density we obtain K = 5.96(4) 
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FIG. 8: Intramolecular distribution function as a function of 



P = r/Rg, for 



.(3) 



(top) and z = z^^' (bottom). We 



report full-monomer (FM-CGR) and tetramer (t) results. 



(tetramer), K = 6.15(3) (FM). For z = z^^) discrepan- 
cies are small (the relative difference is 5% for $ = 6 
and 8% for $ = 10) and so are they for z = z^^) (9% 
for <& = 10). For z = z^'*) and z*^^^ we observe larger dif- 
ferences, the relative difference being larger than 10% 
for <f> > 5 and 3, respectively. Thermodynamic re- 
sults are completely consistent with the structural ones. 
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- FM-CGR, <I>=1 

- FM-CGR+0.2, $=2 

- FM-CGR+0.4, <I>=4 

• CGBM,<I>=1 
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— FM-CGR, $=1 
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FIG. 9: Distribution of i? = Rg,i{$)/(Rg} for z = z'^' (top) 
and z — 2*^^ (bottom). We report full-monomer (FM-CGR) 
and tetramer (t) results. 



FIG. 10: Intermolecular center-of-mass pair distribution func- 



tion as a function of p = r/Rg, for z 



.(3) 



(top) and z = 



.(1) 



(bottom). We report full- monomer (FM), single-blob (SB) 
and tetramer (t) results. 



For z < z^^\ i.e. for systems such that A2 ^ 3 (i.e., 
A2/A2.GS ^ 0.55, where ^2,05 = 5.50 is the good-solvent 
value), the tetramer model reasonably works also deep 
in the semidilute regime. On the other hand, for larger 
values of z, many-blob interactions begin to play a role. 
Hence, in this region, the model is only predictive for 



$ not too large, i.e. for densities such that blob-blob 
overlaps are rare. 

It is also interesting to compare the thermodynamic 
predictions for the SB model. For z = z^^' and z ~ 
z'^) the SB model gives estimates of K that are close 
to the tetramer and FM ones. Differences are observed 
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FIG. 11: Inverse compressibility K as predicted by full- 
monomer (FM) simulations, the single-blob (SB) model (top), 
and the tetramer (t) model (bottom). 



for z > z^^': in this regime the SB model significantly 
underestimates K as long as $ > 1. 



D. Comparison for $ — >^ oo and random-phase 
approximation 

We wish now to compare the thermodynamic behav- 
ior for large values of $. For the polymer system, we 
have K{z,^) « KJ^z^) = /cfmC^)*, where kpAiiz) 
is reported in Sec. IIVBI We now compute the large- $ 
behavior of K{z, $) for the CG models, using the fact 
that the random-phase approximation (RPA)^ is exact 



TABLE IV: Large-$ coefRcients ksB{z), ^4(2), and kpAiiz) 
for five intermediate values of z, see Table[Tl and for the good- 
solvent case (z — 00). For n > 4, k„(z — 00) = 3.204n^~'''^. 



z ksBJz) kjjz) kpMJz) 

z'^> 0.4956 0.501 0.543 

z'-^'> 1.0352 1.083 1.270 

z'^) 1.6271 1.778 2.320 

z'*) 2.2750 2.643 4.072 

z'-^'' 2.9731 3.643 8.504 

00 3.4163 4.584 — 



at large density for systems with soft potentials. For the 
SB model, we start from the virial pressure 






d(3VsBip;z) 
dp 



gcM{p;z,<^)p^dp. 

(28) 

In the RPA, we set gcMif] z, $) = 1, a property which is 
rigorously true for $ — >■ cxd. Integrating by parts we end 
up with the usual RPA expression^ 

ZsB{z,<i>) = l + -<^ksB{z) 



ksB{z)^i l3VsB{p;z)p'dp, (29) 

Jo 

from which we obtain Ksb{z,^) = 1-1- $^5^(2). 
The function ksB{z) is reported in the supplementary 
material^ and in Table ITVl for some specific values of z. 
The RPA expression ([29|) reproduces the SB compress- 
ibility factor at the 1% level for z = z^^^ for all values of 
$. For larger values of z it provides a very good approx- 
imation for $ > 3 (deviations are at most 2% for $ = 3 
and are less than 1% for $ > 8). For smaller densities de- 
viations are larger, due to the fact that the small-density 
behavior is not correctly reproduced (compare ksB{z) 
with the estimates of A2 reported in Table |T| . 

The RPA expression (|29|) can be extended to the 
tetramer model. Starting from the expression of the pres- 
sure in terms of the atomic virial and assuming the ab- 
sence of blob-blob correlations among blobs belonging to 
different tetramers, we obtain for a CGM with n blobs 
(see also Ref. ISTI ) 



Z(z,$) = l+-$A:„(z) 



kn{z) = 3n^ 



/3W{p; z) p^dp. 
(30) 



where W{p; z) is the intermolecular potential. Estimates 
of ^4(2) are reported in Table IIVI In the good-solvent 
case, using the transferability assumption, the integral 
should scale as n^^'^ , so that we obtain kn{z ~ 00) = 
3.204n^~'^'^ (the prefactor has been fixed by using the 
explicit rescaling factor reported in Ref. l30l ) . 

Knowledge of the constants k{z) allows us to compute 
the deviations between CGM and FM results for $ —> 00 



15 



without the need of simulations, since 



Kn{z,<^) _ kn{z) 



K 



'fm{z,^) kpAiiz) 



(31) 



for $ — >■ oo. For z = z^-^\ ksB^z) and fc4(z) are close to 
kFM{z): the SB K{z,^) (or equivalently Z(z,$)) differs 
from the corresponding polymer quantity by 9% for large 
$. For such a small value of z, there is little advantage 
in using the tetramcr model if one is only interested in 
the thermodynamics: the difference decreases only to 8%. 
For z = z'^-* differences are larger: 18% and 15% for the 
SB and tetramer model, respectively. For larger values 
of z, instead, it is quite clear that both the SB and the 
tetramer model grossly underestimate the correct poly- 
mer pressure deep in the semidilute regime. 



V. TRANSFERABILITY IN THE NUMBER OF 
BLOBS 

A. Transferability: general considerations 

In the previous section we have determined a CG 
tetramer model appropriate to describe the 0-to-good- 
solvent crossover. It was found that the predictions of 
the tetramer model are accurate up to a value of the 
density $(z), which increases when z decreases towards 
the 9 region. Above $(z) the tetramer model is observed 
to deviate from the FM results and to approach the RPA 
predictions at large $. We wish now to extend the CG 
model to a larger number n of blobs per chain, which is 
the key ingredient to increase $. In the good-solvent case, 
the basic transferability assumption is that the tetramer 
potentials, expressed in terms of the blob radius of gy- 
ration fg, provide an accurate CGM for any number of 
blobs. For z = oo, the approach is therefore the follow- 
ing. First, we express any tetramer bonding potential 

r/fg-. 



Vij{p; n = 4), p ~ r/Rg, in terms of a 



(32) 



where 7^.4 = fg/Rg for the tetramer. Then, the potentials 
for the n-blob model are defined by 



Vi2{cr;n) = y„_i,„(cr;n) 
%{a-n) = V23{<J-A) 
V,j{a-n) = Vii[a-A) 



^12 (a; 4), 

l*-.7l=2, 

|*-.7|>3. (33) 



Note, that blobs i and j with |j — j| > 3 always inter- 
act with potential Vi4(ct;4) to guarantee the local self- 
repulsion, which is necessary to obtain a good-solvent 
CGM. The angular potentials are instead unchanged 
when increasing the number of blobs per chain. If Rg is 
used as reference length scale, i.e. we consider Vij {p; n) 
with p = r/Rg^ these relations imply the rescalings 
F(p;n) = X^(p7^4/7^„;4)^ 



To justify the previous relations, we should first note 
the dual interpretation of the coarse-graining procedure. 
Up to now, we have considered two blob models with rii 
and 712 blobs each as providing two different represen- 
tations (with different resolutions) of the same polymer 
chain. However, since polymer chains are fractals, hence 
scale invariant, for large degree of polymerization i, the 
multiblob model can be given a different interpretation. 
We assume now that the number m of monomers be- 
longing to a blob is fixed, so that CGMs with ni and 
ri2 blobs are CGRs with the same resolution of polymer 
chains of different lengths, Li = riim and L2 = n2m^ 
respectively. Within this dual interpretation it is easy to 
justify transferability. First, we note that, to a very good 
approximation, the size of the blob depends only on m 
and not on the length of the chain. The latter nontrivial 
property was verified for n > 4 in App. A of Ref. uSk 
Indeed, since Rg ~ L'^ and fg/Rg ^ n~'^ with good pre- 
cision for n > 4, we have fg ~ m'', which depends only 
on m and not on L. Therefore, fg is the same for the two 
CGMs. Second, we assume that the interactions have 
little dependence on the chemical distance between the 
blobs, except for the case in which the blobs are very 
close along the chain, and are insensitive to the length of 
the CGM. If this holds, potentials are the same for the 
two CGMs, i.e. they can be transferred from one model 
to the other one. If we wish to use this result within the 
original CG interpretation, we should simply note that 
potentials are invariant, if fg is the basic length scale. If 
Rg is used instead, a rescaling of the length scale should 
be performed. 

Let us now consider the case of polymer chains in the 
thermal crossover region. Let us consideri, at the same 
temperature, two chemically identical chains of length 
Li and L2, respectively, and their CGRs in terms of ni 
and 712 blobs, each blob consisting of the same num- 
ber 777 of monomers. If we assume that (i) the size of 
the blob is the same for the two chains (this hypoth- 
esis looks reasonable, but we shall show below that it 
is only a rough approximation) and that (ii) interac- 
tions are independent on chemical distance except for 
very close blobs along the chain (again this looks quite 
plausible), we expect the potentials to be approximately 
transferable without any change as in the good-solvent 
case. However, since z = (T — Tg)L^/^, the two chains 
correspond to two different crossover parameters zi and 
Z2- Therefore, the CGM with 77 > 4 blobs obtained by 
using the tetramer potentials at z is a CGM for a system 
at z' > z. In the renormalization-group language, z flows 
towards the stable good-solvent fixed point as the num- 
ber of blobs increases. If hypotheses (i) and (ii) above 
were both approximately correct, we could simply esti- 
mate z' = z(7i/4)-'^/^. However, since z' =^ z, we expect 
fg to differ in the two cases, hence hypothesis (i) docs 
not hold, and therefore the relation between z and z' is 
more complex. 

It is interesting to revisit this argument within the 
usual interpretation of the coarse-graining procedure, in 



16 



TABLE V: The transferability mapping: for each z (2:4) in 
the tetramer model (yl2-4 is the corresponding second-virial 
combination) we report the value of z (zio), the corresponding 
second-virial combination j42-io, and the rescaling factor A 
[see Eq. (|35p ] for the decamer model. 



n = 4 

2:4 



An 



n = 10 

Z\Q ^2-10 



A 



j'^> = 0.056215 0.993 z*^""^' 
z(2' = 0.148726 1.978 2'^""^' 



= 0.08628 1.374 0.9875 
= 0.22953 2.524 0.9833 



2(^' = 0.321650 2.962 z'^''-^*' = 0.47018 3.442 0.9684 
2^ = 0.728877 3.943 2(1''-*) = 1.00862 4.267 0.9604 



.(5) 



2.508280 4.915 2 



(10-5) 



= 3.74117 5.088 0.9778 



which an increase of n corresponds to an increase of 
the resolution of the polymer CGR. In this case z is 
fixed. The previous argument implies that potentials 
should become increasingly softer as the resolution in- 
creases, i.e., the effective z decreases with increasing n. 
This result is completely consistent with the general ar- 
gument of deGennes,— who noted that, close to the Q 
point, polymers show two different spatial regimes. If we 
indicate with Rt the thermal blob size^ blobs such that 
fg <C Rt behave as ideal chains, while excluded-volume 
effects dominate for -fg^Rt. For fg^Rt, which is the 
relevant case here, blobs behave in an intermediate way, 
excluded-volume effects becoming increasingly less rele- 
vant — hence the effective z decreases — as Vg decreases, 
i.e. when the resolution (number of blobs n) increases. 

Given the difficulties presented above, we have devel- 
oped a transferability procedure which works at a formal 
level, without any explicit reference to the underlying 
polymer system. The idea is the following. Consider 
the tetramer set of potentials {V4} at a given value of 
z, z = Z4, and a system of n > 4 blobs interacting with 
these potentials. The question we will ask is: can this 
71-blob system be seen as a CGR of a polymer chain at a 
different value of z, say z„? As we shall show below, the 
answer is positive. We will compute z„ and we will re- 
late the size of the new CGR, i.e. Rg^b^ to the size of the 
corresponding underlying polymer chain. In terms of the 
underlying model, the tetramer potentials appropriate to 
describe a polymer chain at temperature T can be trans- 
ferred without changes to an n blob CGR of the same 
polymer chain, but at a different temperature T' > T. 
In the rcnormalization-group (RG) language, we are con- 
sidering a RG transformation at fixed hare parameters, 
hence we must consider the temperature RG flow towards 
the good-solvent fixed point T = 00. 



B. Definition of tiie higher-resolution models 

Our recipe to transfer the tetramer potentials to ?i-blob 
CGMs works as follows. First, we introduce two universal 
functions: We define FA^b(A2,n) = i32(z)/[i?g,b(z,?T-)]^ 



[i?2 is the second virial coefficient defined in Eq. ([T])], 
which is a universal function of z and n or, equivalently, 
of A2 and n, and 5*6(^42, n) — Rg^h{z,n)/ Rg{z) which is 
also a universal function of A2 and n. Of course, the two 
functions are related by FA.b{A2,'n) = A2Sb{A2,'n)~^. 
These two functions can be computed by FM simulations 
and are assumed known in the procedure we shall present. 
For n = 10, the case we will be interested in, they are 
reported in Appendix [Xj 

The procedure is the following: 

(i) We consider the n-blob model with the same po- 
tentials used for the tetramer. For the bonding 
potentials we set Vij{p;n) = V^b(yo;4), where ij are 
related to ab as in Eq. ([55|. while angular poten- 
tials are unchanged. No rescalings are performed 
at this stage, so that p for the n-blob model should 
not be identified with r/Rg. We will write therefore 
p = r/R, where R simply sets the length scale but is 
otherwise arbitrary (in the numerical computations 
we set R = I). Then, we compute (for instance, by 
Monte Carlo simulations) the second virial coeffi- 
cient B2.MC, the radius of gyration Rg^b.MC, and 
A2,b,MC = B2,MC I Rl,b,MC °f th"^ n-blob chain. 

(ii) We determine the value of A2 appropriate for the n- 
blob model, i.e. A2{zn,n)^ by solving the equation 
FA,b{A2,n) = A2,b,MC- Then, using Eq. (g]), we 
determine the n-blob corresponding z„. 

(iii) The radius of gyration Rg of the polymer chain, 
whose CGR is provided by the n-blob chain, is 
given by Rg^MC = Rg,b,Mc/Sb{A2,n), using the 
value of A2 computed at point (ii) . 

(iv) Once Rg is known, if we wish to express all lengths 
in terms of the radius of gyration, it is enough to 
redefine the bonding and the intcrmolccular poten- 
tials as 



V{p'; z„; n) = V{p'Rg.Mc/R; z; 4), 
where p' = r/Rg{zn; n). 



(34) 



With these definitions, we obtain a model with a higher 
resolution at a different value z„, which has the correct 
A2, hence it gives the correct thermodynamics, and gives 
the correct result for Rg^biz;n)/Rg{z;n), hence it is also 
structurally consistent. 

We have applied this strategy starting from the 
tetramer potentials at z = 2*^^^ . . . , z^'^^. For each of 
them we have determined a decamer model with n = 10 
blobs. The corresponding values of ziq and ^2(210) 
are reported in Table [V] To clarify the procedure, let 
us show how the method works in a specific example, 
applying the transferability procedure to the tetramer 
model at z = z^^> = 0.056215 (Explicit expressions for 
the corresponding potentials are reported in the supple- 
mentary material^^). Consider the decamer model with 
n = 10, using the tetramer potentials (no rescalings are 
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performed). For this model we determine numerically 
[point (i)] B2.McR~^ « 5.740, Rg,b,Mc/R « 1-532, and 
A2,b.MC ~ 1.597. Then [point (ii)], we first solve the 
equation i^^^(,(^2, 10) = 1.597 (i^^_f,(y42, 10) is given in 
App. |X]), obtaining A2 ~ 1.375, and then the equation 
^2(^10) = 1-375 [A2{z) is given in Eq. (g])], obtaining 
^10 ~ 0.0862. For such value of A2 , the results of App. |X] 
give 5fc(A2, 10) « 0.951, hence Rg,Mc/R ~ 1-610. Hence, 
for the dccamcr model, if we set p' = t/ Rg expressing all 
lengths in terms of the radius of gyration of the under- 
lying polymer chain, we should rescale the potentials as 
y(p';10) = F(1.61p';4). In the good-solvent case, the 
rescaling is equal to the ratio TZi/TZn = ^3,4/^9, n- In the 
crossover region, this relation does not hold, because of 
the flow of z. Hence we define 



^^_ rgjZiA) R 

ig,MC 



rgiz4,n) R 



(35) 



which encodes how much of the length rescaling is due to 
the change of the parameter z. At the renormalization- 
group fixed points z = and z = 00, we have 
A = 1. In the example presented above we have 
fg(z4,4)/fg(z4,n) « 1.590, so that A « 0.988. The cor- 
rection is small but not negligible. 

As we discussed in Sec. IV Al a naive application of the 
transferability ideas would predict z„ = 24(71/4)^/^. We 
can check how this approximation works in the present 
case. If we set ziq = Z4(10/4)^/^, for our five values of 
Zi we would obtain zw = 0.0889,0.235,0.506, 1.15,3.97, 
which are close to the estimates reported in Table |Vl 



C. Comparison with full- monomer results 

By definition, our procedure is such that A2 and the 
ratio Rg.b/Rg at zero density are exactly reproduced. We 
wish now to check whether other structural and thermo- 
dynamic properties are satisfied. In Figs. [12] and [13] we 
show the zero-density intramolecular distribution func- 
tion (?intra(/3) and the distribution of the CGR radius of 
gyration i?g,h, respectively. In both cases, the agreement 
is excellent, confirming that the intramolecular structure 
is correctly reproduced. In Fig. [T3] we show the zero- 
density center-of-mass pair distribution function for the 
decamer and for the polymer chain: again, the decamer 
appears to be quite accurate. Note that all these results 
are far from obvious and indicate that our procedure has 
correctly identified the renormalization-group flow from 
the 9 point to the good-solvent fixed point. 

Let us now consider the thermodynamic behavior. By 
construction, the decamer model reproduces the second- 
virial combination A2. Let us now consider the third 
virial combination A3 = B^R^^. In Table FVll we re- 
port the results for the SB model, the tetramer, and 
the decamer model, and compare them with the FM 
predictions ji^ The SB model significantly underestimates 
A3, as already observed in Scc. lIVAl The tetramer model 
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FIG. 12: Intramolecular distribution function gintTa.{p)p^ at 
zero density as a function of p = r/Rg. We report decamer 
(lot) and full-monomer (FM-CGR) results for several values 
of z, see Table JV] For the sake of clarity, results at different 
values of z are shifted upward according to the legend. 
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FIG. 13: Distribution of the ratio R = Rg^j 
density. We report decamer (lOt) and full-monomer (FM- 
CGR) results for several values of 2, see Table IVl For the sake 
of clarity, results at different values of z are shifted upward 
according to the legend. 



appears to be accurate close to the good-solvent limit, 
but some deviations are observed for intermediate val- 
ues of z. The decamer model reproduces the FM results 
with a relative accuracy of less than 1% for all values of z, 
hence 10 blobs are enough to reproduce quite accurately 
the thermodynamic behavior in the low-density regime. 
Finally, let us consider the finite-density behavior. In 
the absence of FM simulations for z'^^^~^^ , . . ., z'-'^°~^\ 
we cannot directly compare structural properties. We 
will thus limit ourselves to compare the thermodynamic 
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FIG. 14: Radial distribution functions for a pair of iso- 
lated chains as a function of the center-of-mass separation 
p — r/Rg. We report decamer (lOt) and full- monomer (FM) 
results at zero density for several values of z, see Table [V] 



TABLE VI: For several values of z, see Table [V] we report 
As for polymers (^Is-fm),— for the decamer (As^io), for the 
tetramer {A3-4), and for the SB model As-sb- 



A3- 



FM 



A3- 



A3- 



A3- 



Z 



j(io-i) 

j{10-2) 
j(10-3) 
^(10-4) 
(10-5) 



0.216 0.217(1) 0.2042(6) 0.152 

1.18 1.191(5) 1.132(4) 0.891 

2.76 2.78(1) 2.540(6) 2.160 

4.95 4.99(2) 4.91(1) 3.944 

7.96 7.95(3) 7.94(1) 6.396 
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TABLE VII: Estimates of A' for the CGM with n = 4, 10 blobs 
and for the single-blob (SB) model. For the polymer model 
we report the field-theory (FT) results of Ref. |5|, see App.lB] 
We also report the relative deviations A„ — 100\Kn/KFT — 1| 
with respect to the field-theory result. Values of z reported 
in Table El 



$ FT n = 10 n = 4 



SB Aio% A4% AsB% 



,{10-5) 



1 4.45 4.42(2) 

2 8.98 8.64(7) 
4 19.56 17.9(2) 
6 31.30 27.5(5) 



4.29(1) 
8.135(25) 
16.18(8) 

23.9(3) 



3.91 0.7 

7.02 3.8 

13.26 8.5 

19.49 12.1 



3.6 12.1 

9.5 21.8 

17.3 32.2 

23.6 37.7 



^(10-3) ^ 2.99 2.965(4) 


2.85(1) 


2.83 


1.0 


4.7 


6.35 


2 5.29 5.18(3) 


4.98(1) 


4.76 


2.1 


5.9 


10.0 


4 10.23 10.00(6) 


9.19(4) 


8.64 


2.2 


10.2 


15.5 


6 15.38 14.6(1) 


13.8(2) 


12.52 


5.0 


10.3 


18.6 



^(10-1) I iQg 1.686(9) 



2 2.41 2.42(2) 
4 3.89 3.91(3) 
6 5.39 5.33(8) 



1.685(4) 
2.386(7) 
3.80(2) 
5.28(7) 



1.67 0.2 

2.36 0.4 

3.75 0.5 

5.15 1.1 



0.3 1.2 

1.0 2 

2.3 3.6 

1.9 4.7 



TABLE VIII: RPA coefficients kssiz) and k„{z), and cor- 
responding full-monomer coefficient kpAiiz); see Sec. IIVDI 
for the definitions. We also report the relative deviations 
Aio — 100\kio {z)/kFM{z) — 1| for the decamer. Values of z 
reported in Table [V] 



z ksBJz) k4{z) fcio(z) kpMJz) Aio(2:)% 



^(lO-l) 
j(10-2) 
j(10-3) 
j(10-4) 
,(10-5) 



0.698 0.727 0.751 

1.357 1.469 1.563 

1.937 2.201 2.397 

2.502 2.976 3.339 

3.103 3.922 4.587 



0.798 


6% 


1.798 


13% 


3.042 


21% 


4.990 


33% 


10.39 


56% 



behavior, using the field-theory expressions^ of App. |Bj 
to compute i^ as a function of $. As we discussed 
in Sec. IIVBI field theory appears to be quite accurate, 
differences from the FM value being at most 4%. In 
Table EIIl we report K{z,^) for the different CGMs 
and compare it with the field-theory prediction. As ex- 
pected, close to the 0-point (z = z^^^^^^), the decamer 
reproduces very precisely the polymer result, at least for 
<& < 6. To estimate the discrepancy for larger values 
of (f> [see Eq. (|5T|) ]. wc use the RPA estimates of the 
asymptotic behavior reported in Table IVIIII Discrepan- 
cies appear to be under control for all values of $, being 
at most 6% in the limit $ -^ 00. For z ~ z^^^"'^^ differ- 
ences are reasonable up to <i> < 6. For larger densities 
the systematic deviations are larger and Z and K are 
underestimated by 21% for $ ^^ 00. For z = z^^^^^^ the 
behavior is similar to that observed in the good-solvent 
regime)^ the decamer model is significantly more precise 
than the tetramer one and appears to be reliable up to 
$«4. 



VI. CONCLUSIONS 

Recently )^^'^" we developed a consistent coarse- 
graining strategy for polymer solutions under good- 
solvent conditions. In this work we extend this strat- 
egy to the thermal crossover region. For large values of 
L, the universal features of the thermal crossover can 
be completely characterized in terms of universal scal- 
ing functions, which can be conveniently computedi^ by 
using the two-parameter model (TPM), at least not too 
close to the 9 point. Taking advantage of this relation, 
full-monomer results have been obtained by performing 
simulations of the Domb- Joyce model;^ which represents 
the lattice version of the TPM. The TPM results are ob- 
tained by varying the on-site repulsion parameter w to- 
gether with the chain length L, in such a way to keep 
the product wL^''^ fixed when taking the scaling limit 
L — > 00. This procedure enables us to explore the ther- 
mal crossover region at fixed z = {T — Tg)L^/^ and re- 
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duccd density <J> and, therefore, provides predictions for 
the scahng functions associated with any generic prop- 
erty of the solution. 

TPM zero-density scaling functions are used as target 
distributions to develop a CG model. As in our previous 
investigation for chains in the good-solvent regime, we 
have first developed a tetramer model, that is a model in 
which each chain is represented by four "blobs." For the 
tetramer model we determine the interaction potentials 
at five different values of the parameter z, that corre- 
spond to different solvent quality. The intramolecular 
potentials are defined in such a way to reproduce the 
structure of an isolated chain. For this purpose we com- 
pute several single-chain scalar structural distributions 
and then use the iterative Boltzmann inversion proce- 
dure to determine the intramolecular potentials. Inter- 
molecular interactions are specified by a single blob-blob 
pair potential which is determined by matching the radial 
distribution function between the centers of mass of the 
chains. The tetramer CG model set up at zero density is 
found to reproduce the collective behavior at finite den- 
sity reasonably well up to a density ^{z) which decreases 
with increasing z. For small values of z, i.e., close to the 
9 regime, the tetramer model can be safely used up to 
very high $ (the error on Z is at most 1% up to $ ~ 10), 
as long as L is large enough to avoid tricritical effects. 
On the other hand, close to the good-solvent regime, we 
recover our previous finding that the tetramer model is 
reliable up to <f> ^ 2. 

In order to enlarge the range of applicability of the 
CG model with density, we have also developed a trans- 
ferability procedure which allows us to use the tetramer 
potentials to build CG models with more blobs per chain. 
While in the good-solvent regime a simple rescaling of the 
characteristic length scale of the tetramer potentials was 
found to be sufficient to provide models with higher CG 
resolutions, in the crossover regime we must both change 
the basic length scale and temperature, i.e., z, as the 
number of blobs is increased. Using this more elaborate 
procedure, we have transferred the tetramer potentials to 
a CG model with 10 blobs per chain. As expected, the 
density range in which the predictions of the decamer 
CG model are accurate is enlarged with respect to the 
tetramer case. 

This work completes our effort to develop a CG strat- 
egy for polymer solutions which employs potentials de- 
rived at zero density, but which is still able to predict the 
correct thermodynamics and structural properties of the 
system at finite density, deep into the semidilute regime. 
Since the potentials are derived at zero density, we avoid 
all inconsistencies related to the use of state-dependent 
potentials;^"— which plague most of the CG models em- 
ployed to study complex fiuids. Our CG model can be 
used to investigate temperature effects in polymer solu- 
tions, hence to compare with experimental results ob- 
tained by using chains of limited extension and that lie 
in this intermediate region of the phase diagram. Exten- 
sions of our strategy to treat polymer-wall and polymer- 



colloid interactions, polymers of different architecture, 
and copolymers are under investigation and will open the 
way to a fully consistent CG modeling of more challeng- 
ing and interesting systems such as polymer solutions of 
various solvent quality in the presence of an absorbing 
wall (depletion interactions), colloid-polymer solutions, 
and block copolymer solutions, the behavior of which is 
well characterized by experiments but not so well by the- 
ory. 



Appendix A: The blob radius of gyration 



TABLE IX: Estimates of the ratios -Rg^tZ-Rg for n = 4 and 
n — 10 blobs for several values of z, see Table [H 



z 


71 = 4 


n = 10 





0.7500 


0.9000 


^(1) 


0.7553(3) 


0.9032(2) 


2*2) 


0.7612(2) 


0.9068(3) 


^(3) 


0.7686(5) 


0.9114(2) 


Z^^^ 


0.7769(5) 


0.9167(2) 


2<^) 


0.7874(4) 


0.9236(1) 


CO 


0.7959(2) 


0.9295(2) 



In this Appendix we wish to compute the function 
Sb{A2,n) = Rg^i,[z,n)/ Rg{z) for n = 4 and n = 10, 
which we parametrize in terms of A2 instead of z. We 
first determine its behavior for A2 — > 0, by performing 
a one-loop computation in the two-parameter model.— >^ 
We begin by considering the average quadratic distance 
between monomers i and j. If L is the total number of 
monomers of the chain, i.e. its contour length, x = ijL 
and y = j/L, we have for j > i: 

^ -((r,-r,)') = 6(y-x) + 



■^3,0 



Hy 



Vlxy 



N3/2 



1/2, 



+ 3(a: 



8:^3/2 (Al) 



_4y3/2 

+4(2 + a; - 3?;)\/l~^ - 8(1 - yfl"^ 



where i?^ q is the radius of gyration for the ideal case 
(2 = 0). We define the swelling factor associated with fg 
as 



Oig(n,z) 






(A2) 



Using the expression reported above we obtain to first 

order in z: 

134 
agil.z) = 1 + Y^^' (A3) 

ag{2,z) = l + -i-(912V2-1181)z 

= 1 + 1.0362, (A4) 

ag{3,z) = 1 + 0.914Z, (A5) 

ag(4,z) = 1 + 0.833Z. (A6) 
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For n > 5, the largc-n approximation 

, , ,'256 1 2 0.98 

ag[n, z) — 1 + z 

works well. Since 



^l,b 



Rl 



We obtain for z ^- 



105 ^Jn n v?!"^ ^ 



Jig 



(A7) 



(A8) 



5-^(^2, 4) = 0.75 + 0.110793Z = 0.75 + 0.00497423^2, 
Sb{A2, 10) = 0.90 + 0.067556Z = 0.90 + 0.00303304^2, 

(A9) 

where we used z — 2^2 (47r)^'^/^ to leading order. To 
obtain the behavior for all values of A2 , we use the results 
reported in Table HXl which have been obtained by FM 
simulations of the DJ model. Interpolating the data we 
obtain 



^6(^2,4) 
^6(^2,10) 



(0.75 + 0.00497423^2 + 0.00023944^^ 



+0.0000679109^i) 



3\l/2 



(AlO) 



(0.9 + 0.00303304^2 + 0.000509387^^ 
-0.000153285^!^ + 0.0000250389^^)^^^ 



Given Sb(n), we can then compute FA,b{n) ~ A2/ Sb(n)" 



Appendix B: Field-theory predictions 



Then, we define (Eqs. (13.27) and (15.12) of Ref.[5|): 



"'(/) = {1~ f)^'-^-^/'"'H{f){l- 0.195 f), (B4) 
2if,w) = /(l-/)-i/(2-)ij„(/)-i/2^HJ^/^^(B5) 



sifiw) 



CqUo 



-(1 - n(3''-2)/wi/ 



uz{f,w) 

^HM)H{fr\l~w^)lw\ (B6) 

where u ^ 8.1075, cq = 1.2, no = 0.53, v = 0.588, 
Lu = 0.80. Here ag is the usual expansion factor, ag = 

Rg/RgO, where Rgo is the value of the radius of gyration 
for the ideal chain, z corresponds to the crossover vari- 
able z apart from a normalization (z = 0.182z), while s 
is cN^^^, where c is the concentration. 

Given z and $, to obtain / and w we work as follows. 
First, we compute / such that z{f, 1) = z/0.182 (/ and 
1/7 = 1 correspond to z and $ = in our variables) and 



3$ 

— ( 

47r 



st = —ag{f) \ 



(B7) 



Then, we determine / and w by solving the equations 



S{f,w) = z/0.182 



s(f,w)=st. (B8) 



The compressibility factor is then determined as (Eqs. 
(17.21) and (17.51) ofRef.i): 



We summarize here the field-theoretical results of 
Schaferf^ reporting the basic formulae which allow one 
to compute the compressibility factor in terms of $ and 
z. In this approach the crossover and density behavior is 
parametrized by two independent variables / and w. The 
variable / parametrizes the crossover from ideal (/ = 0) 
to good-solvent behavior (/ = 1), while w parametrizes 
the density dependence, w = 1 corresponding to <I> = 
and w = to <& = 00. To relate them to z and the poly- 
mer volume fraction $, we need to define several auxiliary 
functions (see Chap. 13 of Ref. Q): 



H{f) - 


= 1-0.005/ -0.028/2 + 0.022/3, 


(Bl) 


Huif) -- 


= (l + 0.824/)°-25^ 


(B2) 


Hnif) = 


= H{f)HuU)-^. 


(B3) 



Z{f.w) = 1 + -(1 



2u*f)WR 



fu*y/^NRQ.. 



\.22Wr 



?,Wr 1 + 1.221^^ 

l + {WR-\){\ + 2WRf/^ 



(B9) 



where u* = 0.364, Wr^ and Nr are functions of w given 

by 



Wr = cono(l - w'^)lw^ 



Nr = no/w\ (BIO) 
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